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ABSTRACT 


NASTRAN analyses were performed with three 
different finite element models on a 90-degree 
pipe elbow to determine principal stresses due to 
internal pressure, inplane bending, out-of-plane 
bending, and, torsion moment loadings. Comparison 
with stresses experimentally obtained under the 
four loading conditions demonstrates the adequacy 
of a finite element model with ideal geometry 
assumptions and an economical mesh spacing. 
Implementation of the NASTRAN modeling is described 
in detail. Deck set-ups for a sequence of computer 
runs to automatically generate NASTRAN input data, 
execute NASTRAN, and postprocess NASTRAN output 
are included. A sample NASTRAN data input file is 
also listed. 


BACKGROUND 


The objective of one of DINSRDC's programs is to reduce the size 
and weight of sea-water piping systems sufficiently to achieve the design 
goals. This objective can be accomplished in several ways: (1) by 
using higher strength materials, (2) by improving fabrication techniques, 
(3) by better designs, and (4) by less conservative structural acceptance 
criteria and stress factors. The third and fourth alternatives require 
the application of more refined structural analysis methods such as the 
finite element method (FEM). With the FEM, the analyst can gain a more 
accurate understanding of general end point motions, which influence 
piping system designs, as well as of motions and stresses in the piping 
systems themselves. 

The design of piping systems depends on a knowledge of the stresses 
and deflections of the pipes due to the anticipated service loads. In 
general, these loads may be either static or dynamic. The static 
category includes pressure, moment, and thermal loads; the dynamic 
includes shock loadings. For any of these loads, the most critical 
locations are usually in components such as elbows and tees. It is well 
known, for example, that a curved pipe does not respond to bending 
loads in the way predicted by elementary beam theory, which is usually 
adequate for straight pipes. 
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This report describes finite element analyses of piping elbows 
subjected to both pressure and moment loads. The NASTRAN finite element 
structural analysis computer program was used. Stresses and flexibility 


factors were computed, and the stresses were compared to experiment. 


The two main purposes of the work were: a 


(1) to determine modeling conditions (including mesh size 
and ideal geometry assumptions) for which a NASTRAN finite 
element analysis of the pipe elbow would yield results in 
good agreement with experiment, and 


(2) to provide deck set-ups for NASTRAN data generation and 


solution implementation which cculd be used to automate 
a computerized parameter study of pipe elbows. 

Such a study, for example, could yield basic information on the 
structural behavior of elbows. Also, given the increased accuracy and 
reliability attainable with finite element models, compared to the less 

\ precise traditional methods, the studies could be conducted in accordance | § 
| with less conservative design rules similar to those of the ASME 


Section LII nuclear pressure vessel code. 


STATEMENT OF THE PROBLEM 


Me Sees eens ow (ieee eee Saree Eee erie 


The structure to be analyzed is the 10-inch, schedule 40, 90-degree 
carbon steel pipe elbow referred to as ME-1 by the Reactor Division of the 


* 
Oak Ridge National Laboratory .+ This elbow, one of four machined at 


a ae Ra S ri SE a Se 


Oak Ridge to test the effects of thinning and ovalling, was machined so 
as to simulate as far as possible an elbow with ideal (uniform) geometry. 


Details on the geometry of ME-l are shown in Figure 1. Oak Ridge also 


*A complete listing of references is given on page 25. 
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made measurements of the actual wall thickness and elbow diameter for 
the test model; these measurements are listed in Tables 1 and 2, 
respectively. All values in each table were averaged to obtain the 
| average wall thickness and outside diameter values of 0.390 inch and 
10.716 inches, respectively. Thus the average value of mean radius (to 
Y middle of wall) is 5.163 inches. 
This elbow structure was subjected to four different loadings: 
(1) internal pressure of 75.53 psi 
(2) in-plane moment of 32,660 in-lb (M, = -32660 in Figure 1) 
(3) out-of-plane moment of 32,660 in-lb (M, = 32660 in Figure 1) 
(4) torsional moment of 32,660 in-1lb (M, = 32660 in Figure 1) 
These loads were also used for the analysis. 
The material properties of carbon steel are E = 2.9 x 107 psi 


(Young's modulus) and v = 0.3 (Poisson's ratio). 


THE NASTRAN MODELS 


Three analyses were performed: 
: Mesh A -~ a uniformly spaced mesh of finite elements (Figure 2) for an 
| ideal (uniform geometry) elbow having the same thickness and 
diameter throughout; the values used for the wall thickness and 
middle radius were the mean values computed from Tables 1 and 2: 
0.390 inci: and 5.163 inches, respectively. 
Mesh B -- similar to Mesh A except that the number of elements in each 
direction was increased by 50%. 
Mesh C -- a uniformly spaced mesh of finite elements for the actual elbow 
ME-1 as measured and tabulated in Tables 1 and 2; the mesh spacing 


was the same as for Mesh A, the coarser of the two idealized models. 


For the ideal models (Meshes A and B), the structure is symmetric 


with respect to the plane Z, = 0 in Figure 1, so that only one-half of the 


structure need be modeled. ne the actual geometry used for Mesh C 
does not exhibit the same symmetry, such symmetry was assumed here also, 
and only half the structure was modeled. This assumption is equivalent to 
assuming that the nonsymmetry in the thickness variation on the unmodeled 


side (25 < 0) has little effect on the solution on the modeled side 
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(2) 20). Since a deliberate attempt was made to machine the ME-1 

model to minimize thickness variations, this is probably a reasonable 
assumption. 
i For all three meshes, the elbow and its adjoining straight sections 
i were modeled with flat plate finite elements (the NASTRAN element QUAD2) 
which have both membrane and bending stiffness. As indicated in Figure l, 
one end was completely fixed, and the other end (where loads are applied) 
was connected to a stiff flange (which was assumed to be rigid in the 
finite element model). 

The bulk of the NASTRAN data deck for each model was generated by a 
computer program called PIPELB, described briefly in the next section. 

The finite element (QUAD2) used in these models yields stresses at 
element centroids. Hence, to facilitate comparison with experiment at the 
middle of the elbow (6 = 45°), an odd number of elements was used in the 
longitudinal direction. 

Figure 2 shows the Mesh A finite element model. This plot was 
generated using the Structural Analysis via Generalized Interactive 
Graphics (STAGING) syateasc 


The characteristics of the three meshes are summarized in Table 3. 


In all three cases, the element aspect ratios within the elbow were 
selected so as not to deviate significantly from unity, which is the 


optimum aspect ratio for the QUAD2 finite element. 


THE DATA GENERATOR 


The computer program PIPELB was written to generate NASTRAN data for ‘ 
THETA-degree elbows of ideal or non-ideal geometry adjoined on both sides 
by straight sections of arbitrary length. In addition to the GRID cards, 
which locate the finite element mesh points, the program also generates 
the plating (NASTRAN data card CQUAD2), pressure loading (PLOAD2), 
symmetry and other constraints (SPC), and flags for congruent elements 
(CNGRNT). As added conveniences, PIPELB will automatically generate the 
set of element identification numbers (ID's) for the mid-THETA elements 
which constitute the solution set for stresses; flexible spokes (CBAR) 


radiating from the centers of both elbow end-circles which facilitate the 


middle of elbow 


Fixed End 


Figure 2 ~ Finite Element Model of Elbow (mesh A) 


SE ES TE TE 


i 
; 
t 
1 


oe 6Etl "39S 6h6 | (0049 309) 2W!2 dd NVULSWN 
61 61 | JUO1JaAeM X14}eW BBes2y | 
| Shlz GHlZ | wopeats yo saaiBep yo saquiny | 
02S 07S | squiod p14Bb yo saquinn | 
| 89% 89h $juaWa|a JO 4aqunu 1eI0L 
| tl il | uo!}2a5 1Yy61e413s | 
15°07 Buruueds Aj[eulpnz!Buo, sqyuawalea jo *oON 
Lt tt | G012DeS 2yH)e3135 
nO°1Z B5uruueds Ajpeutpnz1 Guo, sjzUaWweta yO “ON 
“i Zt Mog [2 a 
906 Buyuueds Aj} eulpnzibuo, sjuawels yo *on 
Zl Zl | uO! IA4IP [elUasaywNII}D Ul SJUaWaLe yO “ON 
| F19e, aes "Ul 06€°0 adid yo ssauyxdsY4t 
Z‘1 Se1qel oas “ul €91°S adid yo snipes ueey 
penzoe jeep! Asqyawoa6 ,en}2e 40 [eap; 


({ 24n6t4y 995) [-3W MOq{2 ediq :pezAleuy Buleg a4njon4ys 


SISATWNY LN3W313 SLINIS 4O AYWWWAS - € 37aVL 


computation of flexibility factors; a rigid element (CRIGD1) at the free 
(loaded) end of the pipe; and the end cap longitudinal load (FORCE) due 


to unit internal pressure. 


The data generator defines grid locations within the elbow section 
in terms of a local cylindrical coordinate system (ID number 1 in Figure 1) 
defined by 
R = RBEND - RPIPE*: cos ¢ 


T= ®, (1) 


Z = RPIPE*: sin 9 


where RBEND is the elbow bend radius measured to the pipe centerline, 
RPIPE is the mean pipe radius, > is the angle spanning the pipe half- 


circumference, and 9, is the angle spanning the elbow arc of THETA degrees. 


1 
(Other coordinate systems used in the NASTRAN data deck are defined in 


Figure 1.) The grid marking system defines the grid ID's as 
GRID ID = 1000 + (500 + ITH) + IPH (2) 


where ITH and IPH are 8, and $, respectively, rounded to the nearest 
integer. The element ID assigned to each quadrilateral is equal to the 
GRID ID that would be assigned for a node located at the element centroid. 

For example, for Mesh A, the elbow was modeled with 17 elements 
covering the 90° elbow arc longitudinally (i.e., a spacing of 5.294°) and 
12 elements covering 180° (i.e., a spacing of 15°) in the circumferential 
direction. Thus, for this mesh, point 51100 is at the intrados (the 
inside edge), 511090 is at the top, and 511180 is at the extrados (the 
outside edge) of the elbow, with all three points at Py = 10.588°. 

The shell was modeled with plate elements (NASTRAN data card CQUAD2). 
Symmetry constraints were applied with SPC's at » = 0° and 4 = 180°. An 
internal pressure loading was specified using a PLOAD2 card for each plate 
element, and the resultant longitudinal load at the closed free end (due to 
the pressure) was specified with a FORCE card. The inplane bending, out- 
of-plane bending, and torsion moment loads were applied at the rigid 
flange at the free end using MOMENT cards. All grid points at the fixed 


end were constrained in all six degrees of freedom. 


Appendix A provides a listing of the card deck needed to create an 
UPDATE program library, compile source code, and catalog executable code 
for the data generator PIPELB. 

Appendix B gives a listing of a typical card deck needed to execute 
the data generator and create a data file suitable for input to NASTRAN. 
The specific example used in Appendix B is for Mesh A. The result of 
executing this step is the NASTRAN data deck, which is listed in Appendix 
E. Various sections of this data deck have been labeled with comment 
cards, which have a $ in card column number one. 

For Mesh C, which uses the measured geometry as listed in Tables 1 
and 2, the data generator interpolates to obtain pipe thicknesses at 
element centroids and to obtain pipe radii at the grid points. This 
interpolation is performed within the generator in subroutine FIND. For 
Mesh C, each element has a unique thickness and hence a unique property 


card (PQUAD2). 


SOLUTION SEQUENCE 


The general computational approach is to perform the following steps: 
Step 1: Compile the elbow data generator, and catalog an executable form 
of the program. See Appendix A. 
Step 2: For the desired mesh, execute the generator, and create and store 
the NASTRAN data deck. See Appendix B. 
Step 3: Using the data deck created in Step 2, execute NASTRAN and save 
(on a permanent file) the stresses at the desired location (8, = 45°). 
See Appendix C for a listing of control cards to accomplish this. 
Step 4: Smooth the stress curves in preparation for plotting, storing 
the results. Figure 3 shows a typical plot of the type desired, 
where various stresses at the section 8) = 45° are plotted as a 
function of 4, the circumferential coordinate. As can be seen in 
Figure 3, the smoothing merely provides a more attractive display of 


the results. The smoothing in Figure 3 was obtained by using the 


B-apline> application subprogram crvrit. “> See Appendix D for a 


listing of the deck to read the NASTRAN principal stresses and per- 


form the smoothing. 
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Step 5: Plot the smoothed results using the informally documented inter- 
active computer graphics package PLOT QUICK developed by Mr. Melvin E. 
Haas of DTNSRDC (Code 1843). 


FLEXIBILITY FACTORS 


The finite element models also allow for the convenient determination 
of flexibility factors, which are of vital interest to piping systems 
designers. To avoid confusion, the term "flexibility factor," as used in 
this report, will be defined precisely. The flexibility factor k of a 
piping component is defined as the ratio of a relative rotation of that 


component to a nominal rotation: 


8 
k= (3) 
where 
a rotation of end "a" of the piping comp nent relative to 
end "b" of that component due to a moment loading M, and 
in the direction of M. 
ae = nominal rotation of an equal length of straight pipe due to 


the moment M. 
For elbows, the nominal rotation is computed using simple beam theory, 


in which case one obtains 


eso (4) 


ae (5) 


for torsional moments, where 
M = applied moment load 
L 


arc length of elbow 
= Young's modulus of material 
shear modulus of material 


= moment of inertia of cross-section 


a me © 
u 


torsional constant of cross-section (equal to the polar moment 
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of inertia for circular cross-sections) 


For 90° elbows, L = 1R/2, where R is the bend radius. 


The determination of the relative rotation 8 ob from the finite 
element model is complicated by the fact that plane cross-sections do not 
remain plane after deformation, so that each point on the surface of the 
shell has a different rotation. However, an “average” rotation for a 
cross-section may be obtained by connecting very flexible beam spokes 
from each shell grid point to a new grid point located at the center of 
the circular cross-section. The rotation of this center point may be used 
as the average rotation for the cross-section. For Mesh A, the center 
point at Station 4 in Figure 1 is labeled 590500 (end "a") and at 
Station 10, 500500 (end "b"). 

Table 4 shows the details of the flexibility factor calculations for 


the three moment loads for Mesh A. 


RESULTS AND DISCUSSION 


Principal stresses were computed by NASTRAN at the middle of the 
elbow (8, = 45°) and compared to experimental results obtained at the 
Oak Ridge National Laboratory. These stresses were computed for four 
separate loading conditions: 
1) internal pressure of 75.53 psi 
2) inplane bending moment of 32,660 in-1b 
3) out-of-plane bending moment of 32,660 in-1lb 
4) torsional moment of 32,660 in-lb 
Figures 4-11 compare the smoothed NASTRAN stress curves as computed 
using Meshes A, B, and C with the experimental data points. The major 
and minor principal stress results are presented on alternating figures 
for each of the four loading conditions. Note that placing more strain 
gauges on the outer surface than on the inner surface of the pipe wall 
resulted in a greater number of experimental data points there. Also, for 
selected experiments, some gauge results were missing. 
The following conclusions may be drawn from Figures 4-11: 
1. The NASTRAN results using Meshes A and B (both ideal geometry) 


are in good overall agreement with the experimental results. Although 
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Figure 4 - Comparison of Computed and Experimental Stresses at 
Middle of Elbow - Major Principal Stress, Pressure Load 
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the coarser model (Mesh A) is non-conservative in that it underestimates 
(in absolute value) the stresses, its prediction of peak principal 
stresses is nevertheless excellent. 

2. The Mesh C (actual geometry) results indicate that, except 
for internal pressure loading, the small variations in geometry shown 
in Tables 2 and 3 play an insignificant role in the stress computations. 
(The Oak Ridge researchers found in their more extensive experimenta! 
that ovality in the pipe cross-section had a major influence on stresses 
only for piessure loading.) 

3. When a large series of computer runs (e.g., a parameter study 
of pipe elbows) is contemplated, the gain in accuracy achieved with 
Mesh B (which has a fifty percent finer mesh spacing than Mesh A) would 
probably not be cost-effective. (The CDC 6400 computer charge at DTNSRDC 
was $292 for Mesh B compared to $99 for Mesh A.) 


ACKNOWLEDGMENTS 


The authors acknowledge with pleasure the fruitful discussions held 
with Mr. L. Kaldor and Dr. Y.P. Lu, both of the Engineering Analysis 
Branch (Code 2744), DINSRDC. 


REFERENCES 


1. Dodge, W.G. and S.E. Bolt, “Experimental Stress Analysis with 
Data Evaluation of Pipe Elbows with Specified Geometric Distortions," 
ONRL-TM-4834, Reactor Division, Oak Ridge National Laboratory, Oak Ridge, 
Tennessee (1978). 


2. Battelle Laboratories, “Interim Report to Air Force Systems 


Command on STructural Analysis via Generalized INteractive Graphics (The 


STAGING System) ," Columbus, Ohio (Aug 1978). 


3. McKee, James M., "B-Spline Functions As a Solution to Some 
Knotty Geometric Modeling Problems," Proceedings of the Symposium on 
Computer Methods in Engineering, University of Southern California 
(Aug 1977). 


4. McKee, J.M. and R.J. Kazden, "G—-Prime B-Spline Manipulation 
Package--Basic Mathematical Subroutines," Report DTNSRDC 77-0036 (Apr 1977). 


5. McKee, James, M., "Updates to the G-Prime B-Spline Manipulation 
Package--26 October 1977," periodic updates available from the author at 
the David W. Taylor Naval Ship Research and Development Center, Bethesda, 
Maryland 20084. 


APPENDIX A 


Listing of Deck to Create UPDATE Program Library, Compile 
Source Code, and Catalog Executable Code for 
Elbow Data Senerator 


(JOB CARD) 

(CHARGE CARD) 

REQUEST ,NEWPL,*PF. 

UPDATE ,F,C.N, 

CATALOG, NEWPL, PFNAME1, ID=XXXX. 

REQUEST,LGO,+*PF. 

FIN, GPT =1,R=3,PL=5000000,1=COMPILE. 

CATALOG, LGO, PFNAME 2, [D=XXXX. 

7/8/9 EOR +*** SGURCE DECK FOR DATA GENERATOR FOLLOWS see 

*DECK ELB 
PROGRAM PIPELB( INPUT .DUTPUT, TAPES2 INPUT, TAPEG=OUTPUT, TAPER. TAPLG, 
~ TAPE10, TAPE11, TAPE12, TAPE13, TAPE14, TAPE15, TAPEI6, TAPE17, TAPE14, 
~ TAPEI9, TAPE20, TAPE21, TAPE22, TAPE23, TAPE24, TAPE25, TAPE26, TAPE27) 


DATA GENERATOR FOR NASTRAN STRESS ANALYSIS OF A PIPE HAVING 
THETA DEGREE ELBOW AND STRAIGHT SECTIONS 
INPUT GENERATION INCLUDES FINITE ELEMENT PLATING (CQUAD2,GRID), 
CONGRUENCY (CNGRNT) AND PRESSURE LOADING (PLOAD2) FOR ELEMENTS) 
AND GRID POINT CONSTRAINTS (SPC) 


MEL MARCUS,DTNSRDC 1843,UAN 1979 


aAagaaAaaaaaa 


DIMENSION IDUM(8), ICNGR( 35,24) 

INTEGER ASET(24) 
COMMON/A/THETA,NPH1,DPH,RPIPS, PIDS, LOADID, PLOAO 
COMMON/B/NX,NY,TABX(9), TABY(7) . TABF(9,7,2) 
INTESER EIO,P1vS,P10,$10,G6(4) 

INTEGER CPC,CDC 

CATA PI/3.141592654/ 


1 FORMAT (1tHT) 
3 FORMAL (10F8.3) 

4 FORMAT (2413) 

8 FORMAT (8A10) 

9 FORMAT (BHCBAR »418,3F8.1,18) 
10 FORMAT (8HCNGRNT ,218) 
11 FORMAT (8HCQUAD2 ,618) 
12 FORMAT (8HGRIO ,218,F8.4,F8.3,F8.4,218) 
15 FORMAT (3HPLOAD2» ,18,F8.3,18) 
16 FORMAT (8HPQUAD2 ,218,F8.3) 
17 FORMAT (8HSPC »318) 
26 FORMAT (8HSET 1 = ,8(16,2H, )/{8X, B(I6,2H, }))} 
27 FORMAT (8X,.16) 
28 FORMAT (SHFORCE ,3t8,F8.4,3F8.0) 


PRESENTED BELOW ARE THE 4 READ ITEMS WHICH OCCUR IN PIPELB 
ALSO INDICATED ARE THE 4 READ ITEMS WHICH FOLLOW 
(IN 2 CALLS TO SUBROUTINE STRATE 


PHI sPORTION OF PIPE CIRCUMFERENCE MODELED, OEGREES 
THETA=PORTION OF ELBOW ARC MODELED ,DEGREES 


AQAaAAANANAAaANAAND 


PHI=160. 
READ (5,3) THETA 
WRITE (6,3) THETA 


NPHI=sNUMBER OF ELEMENTS IN PHI DIRECTION 
NTHETASNUMBER OF ELEMENTS IN THETA OFRECTION (COORD SYSTEM CPC) 


aAaaNna 


READ (5,4) NPHI,NTHETA 


pier anes Wenn Sees 


GAON 


MONANAANAAANAAAAANAAaAN 


ANAAMAAANAAAN 


ann 


40 
45 


WRITE (6,4) NPHI,NTHETA 


RBENO=8END RADIUS MEASURED TO PIPE CENTERLINE 
RPIPS,PIDS=NOMINAL VALUES OF PIPE RADIUS , THICKNESS (PROPERTY ID) 


PI0S=200 
READ (5.3) RBEND,RPIPS 
WRITE (6,3) RBEND,RPIPS 


THIS ITEM CONTAINS A SINGLE READ (BLANK CARD) 
OR MULTIPLE (3+NF*NY) READS 5 
DEPENDING ON WHETHER NF=0 (IDEAL ELBOW) 
OR NF=2 (VARIABLE RADIUS AND THICKNESS) 


NX,NYSNUMBERS OF PHI,THETA-DIRECTION VALUES (TABX, TABY) 
AT WHICH TABULATED VALUES OF 
RADIUS AND THICKNNESS (NF=2) ARE SPECIFIED FOR ELBOW 


IF N&=0, THEN NOMINAL VALUES OF OF PIPE RADIUS 
AND THICKNESS ARE USED THROUGHOUT 
OTHERWISE, VARIABLE RADIUS AND THICKNESS VALUES ARE COMPUTED 
FOR THETA-DIR VALUES BETWEEN TABY(1) ANO TABY(NY) DEGREES 
WITH NOMINAL VALUES USED ELSEWHERE 


READ (5,4) NX,NY,NF 

WRITE (6,4) NX,NY,NF 

IF (NF.EQ.9) GO TO 45 

READ (5,3) (TABX(I),151,NX) 

WRITE (6,3) (TABX(I),1F=1,NX) 

READ (5,3) (TABY(I),1=1,NY) { 
WRITE (6,3) (TABY(I),121,NY) ‘ 
DO 40 K=1,NF 

DO 40 J=1,NY 


READ (5,3) (TABF(I,U,K),1=1,NX) 

WRITE (6,3) (TABF(I,J,K),1=1,NX) 

CONT INUE 

CONTINUE - 
THE 2 ITEMS READ NEXT ARE TRIGGERED BY "CALL STRATE(1)" 
THE ITEMS TO BE READ ARE DESCRIBED ATOP SUBROUTINE STRATE 
THE LAST 2 ITEMS READ ARE TRIGGERED BY "CALL STRATE(2)" 
THE ITEMS TO BE READ ARE DESCRIBED ATOP SUBROUTINE STRATE 
CPC=1 

coc=1 

MID=45 

LOADID=2t 

PLOAD=-1. 


OPtt= PHI /NPHI 
OTH=THETA/NTHETA 
NPH1 =NPHI +1 
NTH1=NTHETA+1 
IHTHeNTH1/2 


END CAP LOAD DUE TO UNIT INTERNAL PRESSURE = 27 


CAPLO=.5*NPHI *RPIPS*RPIPS*SIN(OPH *PI/180. ) 
WRITE (27,28) LOADIO,200500,3,CAPLD,0.,0.,-1. 


28 


GENERATE FE PLATING FOR STRAIGHT SECTION (TOP) 
CHECK SUBROUTINE STRATE FOR INPUT TO BE READ IN 
CALL STRATE (1) 


GENERATE FE PLATING FOR ELBOW 


aNMANA AAaAN 


TH=0. 
IBASE=500000 


on 


DO 200 I=1,NTH1 
PID=PIOS 

RPIPE=RPIPS 

ITH=TH+.5 
INEXT=TH+OTH+.5 
IMOVE=1000*«( INEXT=1TH) 
PH=0. 

1D=I BASE 


DEFINE MYTHICAL CENTER POINTS ALONG ELBOW 


aaa 


NBARS=0 
IDCEN=1D+500 


00 150 J=1,NPHt 


SPOKES (CBAR - 9) CONVERGE TO MYTHICAL CENTERS AT ELBOW ENDS 


a00n aan 


IF (1.NE.1.AND.1.NE.NTH1) GO TO 50 

N@ARS- NBARS+1 

IOBAR= LOCEN+NBARS 

If (J.EQ.1.0R .u.£Q.NPH1) GO TO 52 

WRITE (9,9) 10BAR,PIDS+2,1D0,IDCEN,0.,1.,0.,1 

GO 10 S0 
52 WRITE (9,9)10BAR,PL10S+1,1D, 10CEN,0.,1.,0.,1 
50 CONTINUE 


IPH=PH+.5 
INEXT=PH+O0PH+.5 
JMOVE=INEXT-IPH 


POSSIBLY INTERPOLATE FOR PIPE RADIUS 


aaa 


IF (NF.EQ.0) GO TO 60 
IF (TH.GT.TABY(NY)) GO TO 62 
CALL FIND (PH,TH,1, THICK) 
CALL FINO (PH,TH,2,00) 
RPIPE=.5*(O0-THICK) 
{REG=NTH1 
GO TO 64 

60 IREG=1 
GO TO 64 

62 IREG=I 

64 TREGI=IREG+1 


GRID FOR ELBOW - 13 


AANAD 
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ana 
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aaa 


anNaAANn 


an an 


80 


115 


20 


109 


180 


ELBOW MODELED IN CYLINDRICAL R,7,2 COORD SYSTEM CPC 


IN WHICH T IS PORTION OF THETA 


ANG=PH*PT/180. 
R=RBEND-RPIPE+*COS( ANG) 
Z=RPIPE*SIN( ANG) 

WRITE (13,12) ID,CPC,.R,TH,2,CDC 
IF (J.EQ. NPH1) GO TO 20 

{F (1.EQ. NTHI) GO TO 415 


ICEN=TH+.5+*D0TH+,.5 
JCEN=PH+.5+DPH+.5 
£1D=1000*1CEN+JCEN+500000 
TF (T.EQ.UHTH) ASET(JU)=zEID 


ASSIGNING CONGRUENCY 
ICNGR(1,J)=E10 


POSSIBLY INTERPOLATE FOR PIPE THICKNESS 
1F SO , GENERATE PROPERTY CARD = 16 


IF (NY.EQ.0} GO TQ 80 

IF (1H.GT.TABY(NY)) GO TO 80 
PID=EID 

PC=PH+.S5*DPH 

TC=TH+.5S*OTH 

CALL FIND (PC,TC,1, THICK) 
WRITE (16,16) PIO,MID, THICK 
CONT INUE 


CQUAD2 FOR ELBOW = 11 


G{1)#1D 

G(2)<G(1)+1IMOVE, 

G(3)#G(2)+JMOVE 

G(4)#G(3)-IMOVE 

WRITE (11,11) EID,PID,(G(L),L=1,4) 


PLOAD2 FOR QUADS - 15 


WRITE (15,15) LOADID,PLOAD,EID 
JF (J.NE.1) GO TO 100 


SPC FOR ELBOW - 16,21 
SYM B.C. AT PHi=0,180 - 18 
ANTI= SYM B.C. AT PHI=0,180 — at 


SI0=71 
L=345 
WRITE (18,17) SID,IO,L 
S10=72 
L=126 
WRITE (21,17) SID,ID,t 


1D=1O+UMOVE 
PH=PH+DPH 


CONTINUE 


TBASEXIBASE+ IMOVE 


30 


RK eres yen 


TH =TH +0TH 
200 CONTINUE 


CNGRNT (THETA DIRECTION ONLY) FOR ELBOW = 10 


aAaAan 


IF (EREG .GE.NTHETA) GO TO 210 

00 205 J=1,NPHI 

DO 205 I=IREG!,NTHETA 
205 wRITE (10,10) ICNGR(IREG, JU), ICNGR(I,u) 
2tO CONTINUE 


SOLUTION SET FOR CASE CONTROL ~ 26 


aaa 


WRITE (26,26) (ASET(JU),J=1,NPHI) 
WRITE (26,27) 999999 


GENERATE FE PLATING FOR STRAIGHT SECTION (BOTTOM) 
CHECK SUBROUTINE STRATE FOR INPUT TO BE READ IN 
CALL STRATE (2) 


AO ANNnN 


WRITE (6,1) 
DO 1000 1=9,27 
REWIND [ 

919 READ (1,8) (IDUM(K),K=1,8) 
IF (€0F(I).NE.0) GO TO 920 
WRITE(6,8) ( 1OUM(K) ,K=1,8) 
GO TO 910 

920 REWIND I 

1000 CONTINUE 


; REWLNOD 98 
Cc 
c 
4 sroe 
END 
; SUBROUTINE STRATE (IFLAG) 
4h Cc 
i Cc GENERATE FE PLATING FOR STRAIGHT SECTION 
c IFLAG=1,2 INDICATES PORTION ABOVE,BELOW ELBOW ! 
Cc 3 
OIMENSION 2S( 25), 1Z0Z2CH(6} ,KOUNT( 5), TUCNGR(S, 500) 
COMMON/A/THETA,NPH1, DPH, RPIPS, PIDS, LOADID, PLOAD 
INTEGER CPS,CDS,EID,PIDS.PID,S10,G(4) 
Cc 
2 FORMAT (18,9F8.2/(10F8.2)) 
4 FORMAT (2413) 
10 FORMAT (8HCNGRNT ,218) 
11 FORMAT (8HCQUAD2~ ,618) 
12 FORMAT (8HGRID ,218,F8.4,F8.3,F8.4,218) 
15 FORMAT (8HPLOAD2 ,18,F8.3,I8) 
17 FORMAT (8HSPC ,318) 
25 FORMAT (8HCRIGDI ,318) 
c 
Cc 
c STRAIGHT SECTION MODELED IN CYLINDRICAL R,T,2Z COORD SYSTEM CPS 
Cc IN WHICH Z 1S DISTANCE ALONG STRAIGHT EDGE 
c ( Z=0 AT POINTS OF CURVATURE/TANGENCY ) 
c 
c 
c 


pita ste 


31 


eaana AaAANAA AHAAAMAAN 


an 


aan 


650 


20 


30 


N2SteNUMBER OF GRIO POINTS IN 2 DIRECTION 
25(1)=Z-COOROS OF THE NZS1 PTS IN COORD SYSTEM cPS 


FOR IFLAG=1, VALUES BEGIN AT Z*-Lt AND END AT 220 
FOR IFLAG=2, VALUES GEGIN AT Z20 AND END AT 2212 
WHERE L1,L2 ARE PIPE LENGTHS ABOVE, BELOW ELBOW 


READ (5,2) NZS14,(2S(1),1=1,NZ$1) 
WRITE (6.2) NZS1,(ZS(1),121,NZS1) 


NOZCHS=NUMBER OF CHANGES IN DZ2ZS(141)-zS(1) | 
TZDZCH(K)SVALUES OF I ,I=1 TO NZS1, AT WHICH DZ CHANGES 
(CONVENTION IS TO SPECIFY 121 AS A CHANGE) 


READ (5,4) NOZCHS, (1ZDZCH(K),K#1,NDZCHS) 
WRITE (6,4) NOZCHS, (1ZDZCH(K) , Ket ,NOZCHS) 


SET UP CONGRUENCY SETS 


1B00M=0 

NSETS= NDZCHS 
00 650 I=1,NSETS 
KOUNT(1I)=0 


IF (IFLAG.EQ.2) GO TO 20 
CPS=3 

cOS=3 

NTPGRO=12 

NTPSYM=17 

NTPAS=20 

ILIM=NZS1-1 
IBASC=(500-10*ILIM)*1000 
IRIGIO=1BASE 


GO TO 30 
CPS=2 

cOS=2 

NTPGRD=14 
NEPSYM=19 
NTPAS=22 
TLIM:NZS1 


IBASE=(S00+THETA) #1000 
RPIPE=RPIPS 
PIO=PIOS 


DO 300 I=1,ILIM 
IMOVE=10000 
PH=0, 

IO0=I BASE 


CO 350 J=1,NPHI 
TPH=PH+.5 

INEXT=PH+DPH+.5 
JMOVE=INEXT-IPH 


RIGID ELEMENTS AT FREE END ~ 25 


IF (1.NE.1.OR.IFLAG.NE.1) GO TO 27 
WRITE (25,25) JU+#1,200500, IRIGIO 


aaa aaaa 


aaa 


aan 
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27 


202 


700 


204 


215 


220 


IRIGID=IRIGID+JUMOVE 
CONT INUE 


GRID FOR STRAIGHT SECTION - NTPGRO 


IF (IFLAG.EQ.2.AND.1.EQ.1) GO TO 202 
NT=NIPGRO 
WRITE (NT,12) 10,CPS,RPIPE,PH, 278(1),COS 


TF (J.EQ. NPH1) GO TO 220 
IF (IFLAG.EQ.2.AND.1.EQ.I1LIM) GO TO 215 
EID=ID 


IF (IFLAG.EQ.1.AND.1.E€Q.ILIM) GO TO 204 
IF (TFLAG.EQ.2.AND.1.EQ.1) GO TO 204 


TESTING FOR AND ASSIGNING CONGRUENCY 


1 BOOM=1B00M+41 

I Z2D2CH{NDZCHS+1 ) =NZS1 

DO 700 K=1,NOZCHS 

NSET=0 

1Z1=120ZCH(K )-1 

17221 70Z2CH(K41)=-1 

IF (1.GT.1Z21.AND.1.LE.1Z2) NSET#K 
IF (NSET.EQ.0) GO TO 700 
{F (1800M.EQ.1) NSETIENSET 
NSETS=NSET 
KOUNT(NSET)=KOUNT(NSET) +1 
KNT=KOUNT (NSET ) 
IUCNGR(NSET,KNT) 2EID 
CONTINUE 


CQUAD2 FOR STRAIGHT SECTION = 11 


G(1)=10 

G(2)=G(1)+IMOVE 

G(3)=G(2)+UMOVE 

G(4)=G(3)-IMOVE 

WRITE (14,11) EIO,PID,(G(L),L=1,4) 


PLOAD? FOR QUADS - 15 


WRITE 415,15) LOADICD,PLOAD,EID 

IF (IF LAG.EQ.1,AND.1.€Q.1) GO TO 300 

IF (uU.EQ.1)} GO TO 220 

IF (ITFLAG.EQ.2.AND.1.EQ.I1LIM) GO TO 225 
co TO 300 


SPC FOR STRAIGHT SECTION - NTPSYM,NTPAS, 23,24 
SYM B.C. AT PHI=0,180 - NTPSYM 

ANTI~- SYM 8.C. AT PHI=0,180 - NTPAS 

ELIM ZERO DOF AT PHI=0,180 - 23 

CONSTRAIN FIXED END - 24 


IF (1.€Q.1) GO TO 300 

S$10=71 

L=246 . 
NT=NTPSYM 

WRITE (NT,17) SIO,IO,L 

S$1D=72 

L=135 

NT=NTPAS 


225 


oan 


350 


400 


anon Aan 


405 


an 


agANAN 


10 


15 
20 


30 


35 


40 


50 


WRITE (NT,17) SIO,ID,L 

$10=80 

is4 

WRITE (23,17) SID,L0,L 

TF (1FLAG.EQ.2.AND.1.EQ.1LIM) GO TO 225 
GO tO 300 

S$10=85 

L=123456 

WRITE (24,17) SI0,10,4 


10=1D+UMOVE 
PH=PH+DPH 


CONTINUE 
IBASE=1BASE+IMOVE 


CONT INUE 


CNGRNT (BOTH OLRECTIONS) FOR STRAIGHT SECTION - 10 


DO 405 K=NSET1,NSETS 

KNTSKOUNT(K) 

DO 405 L=2,KNT 

WRITE (10,10) IUCNGR(K,1),IUCNGR(K,L) 


RETURN 
END 
SUBROUTINE FIND (X,Y,K,FXY) 


INTERPOLATE FOR FUNCTIONAL VALUES WHICH ARE 
TABULATED FOR BOTH XK AND Y TABULATED VALUES 


K INOICATES FUNCTION BEING CONSIOERED 


OIMENSION FX(7) 
COMMON/B/NX,NY,TABX(9), TABY( 7), TABF(9,7,2) 


DO 10 J=2,NY 

{F (TABY(U).GT.¥) GO TO 15 
CONTINUE 

U2=NY 

GO TO 20 

J2=eJd 

Jtzl2~1 


CO 30 1=2?,NX 

IF (TABX(1).GT.X) GO TO 35 
CONTINUE 

I12=NX 

SO TO 40 

T2=1 

Tt=f2-1 


XRAT=(KX~TABX(11))/( TABX(12)~TABA(11)) 

00 50 JzJ1,u2 

FRCS) TABE( 11,5, K)tXRAT©( TABF(I2,U,Kj}-TABE(T1,U,%)) 
YRAT#(Y¥~-TABY(J1))/( TABY(JU2)-TABY(J1)) 

FRYZEXC U1) +VRAT# CFX( UZ)“ FR( U1) ) 


RETURN 
END 


6/7/8/3 EOF 
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APPENDIX B 


Listing of Deck to Execute Elbow Data Generator 
and Create NASTRAN Input Data File for Mesh A 


: (uOB CARD) 

3 (CHARGE CARD) 

ATTACH, PIPELB, PFNAME2, ID=XXXX. 

PIPELB. 

REQUEST ,NEWPL. «PF. 

UPDATE,F,C.N,D. 

CATALOG, NEWPL, PFNAME3, ID=XXXX. 

COPYSBF,COMPILE. 

7/8/9 EOR *** INPUT DATA (MESH A) FOR DATA GENERATOR FOLLOWS #206 


90. 
12 17 é 
15.9 5.163 
(BLANK CARD) 
12 -21. -18.5 “16.5 -14.5 “12.5 -10.5 “8.5 “6.5 
-2.5 “1. QO. 
& 1 2 10 11 
12 0. 1. 2.5 4.5 6.5 8.5 10.5 12.5 
16.5 16.5 20.5 
3 1 2 3 
7/8/93 EOR x*ae INPUT DECK TO UPDATE FOLLOWS #48 
*DECK DATA 


NASTRAN CONF1G=6,FILES2(PLT2,NPTP,OPTP, INPT) 
ID PIPE,ELBOW 
APP DISP 
SOL 1,0 
TIME 25 
SSEQUENCE YES 
$GRID=525 cides 
$CONFIG=6 
CEND 
TITLE =STRESS ANALYSIS OF ELBOW ME-t (ORNL~TM-4834) 
MAXLINES=500000 
*READ FAPE26 
DISP=ALL 
STRESS( PRINT, PUNCH) =1 
SPCFORCE=ALL 
$ 
$ SUBCASES HAVING SAME B.C. ARE GROUPED TOGETHER 
$ TO AVOID UNNECESSARY MATRIX OECOMPOSITIONS 
$ 
SUBCASE 1 
SUBTITLESUNIT PRESSURE LOAD 
LABEL=SYM B.C. 
LOAD=21 
SsPCc-91 
OLOAD=ALL 
STRESS=1 
SUBCASE 2 
SUBTITLESUNIT IN-PLANE (2) MOMENT 
LABEL=SYM B.C. 
LOAD=22 
SPC=91 
OLOAD=ALL 
STRESS=1 
SUBCASE 3 
SUBTITLE=UNIT QUT-OF=PLANE (X) MOMENT 
LABEL=ANTI-SyYm B.C. 
LOAD=23 
$PC=92 
OLOAD=ALL 
STRESS=1 
SUBCASE 4 
SUBTITLEZUNIT TORSIONAL (¥) MOMENT 


LABLL=ANTI-SYM 8.C. 

! LOAD=24 

SPC=92 
OLOAD=ALL 
STRESS=1 

SUBCOM 21 
SUBTITLE=PRESSURE LOAD OF 75.53 PSI 
SUBSEQ=75.53,0.,0.,0. 

SUBCOM 22 
SUBTITLE=IN-PLANE MOMENT OF 32660 IN-LB 
SUBSEQ=0. ,-32660.,0.,0. 

SUBCOM 23 
SUBTITLE=OUT-OF-PLANE MOMENT OF 32660 IN-LB 
SUBSEQ=0.,0.,32660.,0. 

SUBCOM 24 
SUBTITLE= TORSIONAL MOMENT OF 32660 IN-LB 
SUBSEQ=0.,0.,0. , 32660. 


STRUCTURE ». PIPE (RPIPE=5.163) CONSISTS OF 
90 DEGREE ELBOW (RBEND=15.0) 
ADJOINED BY STRAIGHT SECTIONS 


SYMMETRY «IMPOSED AT PIPE HALF-CIRCUM (PHI=0,180) 


COORD SYSTEMS...(ELBOW) R=RBEND-RPIPE*COS( PHI) 
T=THETA, (THETA=0,90 AT P.C.,P.T.) 
Z=RPIPE*SIN( PHI) 


COORD SYSTEMS...(STRAIGHT SEC) R=RPIPE 
T=PHI 
(TOP; Z=-DISTANCE FROM P.C. 
(BOTTOM) Z=DISTANCE FROM P.T. 


NODAL VALUES ... (ELBOW) 1D=1000*(500+1TH)+I PH 
WHERE ITH ARE ROUNDED NODAL VALUES OF THETA 
IPH ARE ROUNDED NODAL VALUES OF PHI 


(STRAIGHT SEC) 
(TOP) 10#1000*(500-10*(JU-I))+1PH 
120,1,2,...,U-1 
(BOTTOM! [021000*(5904+101 )+IPH 
T=1,2,..+,.K 
WHERE J IS NO. OF Z-INTERVALS SUBDIVIDING TOP SEC 
K IS NO. OF Z-INTERVALS SUBOIVIDING BOT SEC 


EGIN BULK 
SPOKES AT ELBOW ENDS CONVERGE TO MYTHICAL CENTER POINTS 
READ TAPES 


ELEMENT CONGRUENCY...( ELBOW) THETA DIRECTION ONLY 
(STRAIGHT SEC) BOTH DIRECTIONS 


READ TAPE1O 


ELBOW OEFINED IN COORD SYSTEM 1 
STRAIGHT BOTTOM DEFINED IN COORD SYSTEM 2 
STRAIGHT FOP DEFINED IN COORD SYSTEM 3 
ORNL CARTESIAN SYSTEM IS NO. 4 


PPAABP EAHAAR =~ HHA TD HANG AHA AK AH HAA GH AAA AAA AHAHAHAAMMHAHAAAHAMH 
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$ 


CorRD2Cc 1 t°) 0. QO. Oo. 0. 0. 1. +COR! 
+CORI 1. 0. oO. 

CORD2C 2 t¢] QO. 15. 0. -1, 15. 0. +CoR2 
+COR2 0. 14. 0. 

corp2c 3 0 18. 0. 0. 15. 1. oO. +COR3 
+COR3 14, QO. 0. 

CORD2R 4 Q QO. Qo. 0. 0. 0. “1. +COR4 
+COR4 1. Q. Q. 

$ 

$ FINITE ELEMENT PLATING 

$ 

*READ TAPE11 

$ 

$ RIGID CONNECTION FOR FREE (.OADED) END OF PIPE 

$ 

CRIGDI 1 200500 100500 


*READ TAPE25 
END CAP LOAD DUE TO UNIT INTERNAL PRESSURE 


READ TAPE27 


GRID POINTS 


MYTHICAL CENTER POINT 1 UNIT BEYOND FREE END OF PIPE 


ps) 
~ 
o 


100500 Oo 15.0 -22.0 Qo. 4 


MYTHICAL CENTER POINT AT FREE END OF PIPE 


az 
—_ 
o 


200500 0 15.0 -21.0 0. 4 
STRAIGHT SECTION ABOVE THETAsO 
READ TAPE12 , 


MYTHICAL CENTER POINT AT THETA20 END OF ELBOW 
(EQUIV TO ORIGIN OF COORD SYS 3) 


és) 
i 


0 500500 0 15.0 0. 0. 4 
ELBOW POINTS 
READ TAPE13 
MYTHICAL CENTER POINT AT THETA#90 END OF ELBOW 
(EQUIV TO ORIGIN OF COORD SyY§S 2) 
RID 590500 0 0. 15. 0. 4 
STRAIGHT SECTION BELOW THETA=90 


READ TAPE14 


> 
a) 
_ 


45 2.947 0.3 7.324-4 STEEL 


BENDING MOMENTS APPLIED AT FREE ENO 


PAA ZA *tAAGAO HAHA HAHA OHAMKH tt HRHHOHHUEOKHHHAMHK*AHOM 


MOMENT 22 200500 0 0.50 0. QO. 1. 
MOMENT 23 200500 0 0.50 t. 0. QO. 
MOMENT 24 200500 © 0.50 0. 1. Q. 
$ 

PBAR 201 45 1.0-6 1.0-9 1.0-9 2.0-9 
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PBAR 202 45 2.0-6 2.0-9 2.0-9 4.0-9 


$ 

$ INTERNAL PRESSURE LOADING FOR PIPE j 
$ 
*READ TAPE1S 
{ $ . : 
: PQUAD2 200 45 0.390 i 
: $ : 
$ SPC...(S1D=71) SYM B.C. AT MIDPLANE CUTTING PIPE CIRCUM : 

$ ; 
§ SPC 71 200500 345 : F 

*#READ TAPE17 : ! 

spc 71 500500 345 i 
i *READ TAPE19 ‘ 
SPC 71 590500 345 : 
q *READ TAPE19 : 
: $ 

$s (SI0=72) ANTI-SYM B.C. AT MIDPLANE CUTTING PIPE CIRCUM 

$ 
: sec 72 200500 126 
: *READ TAPE20 

SPC 72 §00500 126 

*READ TAPE2t 

SPC 72 590500 126 

*REAO TAPE22 

* $ 

$ ELIMINATION OF ZERO STIFFNESS DOF 

$ 

*READ TAPE23 

7 i 

$ CONSTRAIN FIXED END 

$ 

READ TAPE24 

s 

Shca0D 91 "1 80 85 

2PCA0D 92 72 80 8S 

ENODATA : 

6/7/8.9 EQF : : 
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APPENDIX C 


: Listing of Deck to Execute NASTRAN and 
Save Stresses at Middle of Elbow 


(JOB CARD) 

(CHARGE CARD) 

f LIMIT, 13000. 

MAP, OFF, 

ATTACH, OLDPL, PFNAME3, 1D=XXXX. 
UPDATE ,F ,C=DATA,N,D. 

ATTACH ,NASTRAN. 

REQUEST ,PUN,*PF. 

RFL, 160000. 

BEGIN, NASTRAN, NASTRAN, , DATA, , PUN. 
CATALOG, PUN, PFNAMEG, I D=XXXX. 


7/8/9 EOR  **** UPDATE CORRECTIONS (IF ANY) FOLLOW **** 
6/7/8/9 EOF 


ee 


APPENDIX D 


Listing of Deck to Smooth Principal Stresses and 
Create Solution File in PLOT QUICK Format 


(JOB CARD) 
‘CHARGE CARD) 
ATTACH, FAPE9, PFNAME4S, ID=XXXX. 
REQUEST, TAPEtO, *PF. 
ATTACH, PROD, 461PRODUCTS, IDacSyS. 
ATTACH, BSPLNLB, 1D=CAMK. 
LIBRARY, PROD, BSPLNLB. 
LOADLOR. 
FIN. 
LOADLOR. 
tGOo. P 
CATALOG, TAPE10, PFNAMES, ID=XXXX. 
7/8/93 EOR #*ee SOURCE DECK FOLLOWS **#8 
PROGRAM PREPLT( INPUT, OUTPUT, TAPES= INPUT, TAPEG=OUTPUT, TAPE, TAPE10) 


THIS PROGRAM READS (TAPE9) NASTRAN-GENERATED ELBOW STRESS OuT®uT, 
SELECTS FROM IT A SOLUTION SUBSET, A 
SMOOTHS THIS SOLUTION SET, 

AND THEN SCALES + SAVES THE ORIGINAL AND SMOOTHED SOLUTION SETS, 
EACH IN A PLOTQUICK (INTERACTIVE PLOTTING) FORMAT ON TAPE10 


CARD INPUT CONSISTS GF NPTS,THE NO. OF ELMTS IN CIRCUMFER DIREC, 
FOLLOWEO BY A SHORT AND LONG NAME, TITS AND TITL, RESPECTIVELY, 
FOR EACH OF THE 34 CURVES WHICH GIVE, 

IN AN ALTERNATING (UNSMGOTHED AND SMOOTHED) PATTERN, 

THE CIRCUMFERENTIAL COORDINATES, FOLLOWED By A CONSECUTIVE STRING 
OF MAJOR-OUTER,MINOR-OQUTER,MAJOR@-INNER,MINOR-INNER STRESS VALUES 
FOR EACH OF THE 4 DIFFERENT LOADING CONDITIONS 


MEL MARCUS, DTNSROC 1643, FEB 1979 


ANNAANANNANNANNAANAANNAAANANAA 


DIMENSION 2(19),TITS(1),TITL(S) 
DIMENSION GPVAL(16),.X(18)5 NSEEK(4),C1(2,18) ,C2(2,103),¥(4,4,18) 
COMMON wORK( 3000) 


FORMAT(1H1) 

FORMAT (13F6. 2) 

FORMAT(12G11.4) 

FORMAT (8A10) 

FORMAT(1X,8A10) 

FORMAT (2613) 
FORMAT(110,8X,3E18.6/(18X,3E18.6,8X) ) 
FORMAT (213,€16.8) 


WwOrnYNaWA — 


NLOADS=4 
NSKIP=7 
NVALUS=16 
NS=4 
NSEEK(1)=6 
NSEEK(2)=7 
NSEEK(3)2=14 
NSEEK(4)215 


READ (5,9) NPTS 
WRITE (G,9) NPTS 
DANG=180./NPTS 
N1#6+(NPTS—1 )+1 


ins ipl nc 


Serpe az 


WRITE (6,1) 
DO 20 K=t,NLOAOS 
WRITE (6.9) K 


DO 10 I2t,NSKIP 
10 READ (9,4) TITS(1) 


DO 20 I=t,NPTS 
REAO (9.8) 1D, (GPVAL(M) ,Mz1,NVALUS) - 
j WRITE (6,8) ID, (GPVAL(M) ,M=1,NVALUS) 
: X(1)s( 17.5) *OANG 


] 
i 00 20 J=1,NS 


N=NSEEK(U) 
20 Y(K,J,1)=GPVAL(N) 


SCFAC=.001 


a aan 


WRITE (6,1) 
KOUNT=90 

DO 100 K=1,NLOADS 
DO 100 J=t,NS 


po 30 I=1,NPTS =, 
C1(1,1)=X(1) 
30 €1(2,1)=¥(K, ue 1) 


CALL CRVFIT(C1,2,NPTS,1,1-0,3,0+1-E"6,N1 0, WORK, 3000 ,C2,RMS, TFAIL) 


IF (K.NE.1.0R.JU.NE.1) GO TO 40 
; KOUNT=KOUNT+1 
7 WRITE (6,7) KOUNT 
READ (5,4) TITS(1) 
wRITE (6,5) TITS(1) 
WRITE (10) TITS(1) 
READ (5,4) (TITL(1),121,8) 
WRITE (6,5) (TITL(I), 11,8) 
WRITE (10) (TITL(I),1=1,8) 
WRITE (6,3) (C1(1,1),121,NPTS) 
WRITE (10) (C1(9,1),129,NPTS) 
KOUNT=KOUNT+ 3 
WRITE (6,7) KOUNT 
READ (5,4) TITS(1) 
WRITE (6,5) TITS(1) 
WRITE (10) TITS(1) 
READ (5,4) (TITL(I),121,8) 
WRITE (6,5) (TITL(1), 121,98) 
WRITE (10) (TITL(T), 11,8) 
WRITE (6.3) (C2(1,1),181,N1) 
: WRITE (10) (C2(1,1),181,N1) 
40 CONTINUE 


KOUNT=KOUNT+1 
WRITE (6,7) KOUNT ~ 
READ (5,4) TITS(1) 
5) TITS(1) 
) TITS(1) 
4) (TITU(E), 121.8) 
5) (TITU(T),181,8) 


( 
WRITE (1 
READ (5 
WRITE (6 


g 
! WRITE 


ear 


WRITE (10) (TITLE), 221.8) 
WRITE (6,3) (C1(2,1)*SCFAC,1=1,NPTS) 
WRITE (10) (C1(2,1)*SCFAC,121,NPTS) 
KOUNT=KOUNT+ 1 

WRITE (6,9) KOUNT, 1FAIL,RMS 

READ (5,4) TITS(1) 

WRITE (6,5) TITS(1) 

WRITE (10) TITS(1) 

READ (5,4) (TITL(1),1=1,8) 

WRITE (6,5) (TITL(1L),121,8) 

WRITE (10) (TITL(I),121,8) 

WRITE (6,3) (C2(2,1)*SCFAC,1=t,N1) 
WRITE (10) (C2(2,1)*SCFAC,1=1,N1) 


c 
100 CONTINUE 
c 
c 
c 
END 
7/8/9 EOR wees INPUT DATA (MESH A) FOLLOWS eens 
12 
PHI--U 
PLPE CIRCUMFERENCE (OE€GREES) 
PHI 
PIPE CIRCUMFERENCE (DEGREES) 
PRES-UOMUR 
NASTRAN MAJOR STRESS AT OUTER STATION 7 OUE TO PRESSURE LOAD 
PRES-NOMUR 


NASTRAN MAJOR STRESS AT OUTER STATION 7 OUE TO PRESSURE LOAD 
seeeessenee SHORT AND LONG NAMES FOR 30 MORE CURVES FOLLOW *etenucereens 
6/7/6/9 EOF 
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APPENDIX E 


Listing of NASTRAN Input Data File 
for Mesh A 


crs 


viva 
vivo 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
vivo 
viva 
vivo 
vivd 
yiva 
viva 
viva 
viva 
viva 
vivd 
viva 
viva 
viva 
vivd 
viva 
vivd 
viva 
vivgd 
viva 
vivo 
viva 
viva 
vivd 
vivo 
viva 
yvavd 
viva 
viva 
vivd 
viva 
vivd 
vivd 
viva 
viva 
vivd 
vivd 
vivo 


zz WOOSNS 

"0° °0'° CES’ SL=03SENS 

ISd €S°SL 40 QvOl 3YuNSS3ud=37L1 LENS 
1% WODENS 

b=SS3M1S 

Viv=ayv010 

Z76=9dS 

vZ=Gv0r 

°O°S WAS-ITANV=138V1 

INSWOW (A) TWNOISYOL LINN=3TLII€8NS 
v 3SvoENs 

b=SS341S 

J1w=0v010 

76=9dS 

€Z=0vo01 

°3°@ WAS-TAINV=T39"7 

AN3WOW (X) 3NWId—40-1NO LINN=3TL115NS 
€ 3Svo8Ens 

b=SS3mM1S 

V1¥=0v010 

b629dS 

ZZ=0v01 

"3° WAS=129V1 

AN3WOW (Z) 3NVId-NI LINN=37151LeNS 
@ asvoens 

t=SS341S 

T1v=0vC10 

$6=9dS 

bZ=QVC1 

°O°S WAS=1398VT 

QvO1 3YuNSS3Yd LINN=37L118NS 
tb 3aSvoENS 
$s 
SNOILISOdNOD3O XIYLYW AYNVYSSZD3NNN OIOAY O1 $ 
¥3H13901 G3dNOYD Juv ~9°a 3IWVS ONIAVH S3SVIENS $ 
$s 
VIW=39a059dS 
b=(HONNc ANTUd )SS3usS 
Viv=dS10 

666666 
“ELISPS ‘8SISPS ‘EpiSeS ‘EcISrs 


“ELISDS ‘S860SPS ‘E80SPS ‘B90GHS ‘ESOSHS ‘8E0SHS ‘EZOSGHS ‘e00SPS = i 13S 


00000S=S3NITKYW 

(veSb-wi-1NYUO) 1L-3W MOB13 40 SISATYNY SS3¥iS= 3ATLIL 
GN39 

9=914NOD$ 

S2S20180$ 

S3A 3ONINOISS 

SZ 3wll 

os 10S 

dSI0 dav 

MOG13‘'3aId OI 

(LONI SdldO’dldN 211d) £S3113‘°9=9I5N0D NWYLSYN 


Ses ec es . a ig See a ai es SS ee eee — 
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viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
viva 
yivad 
viva 
viva 
viva 
viva 
viva 
vivd 
vivd 
yvaivd 
viva 
viva 
viva 
vivd 
yawd 
vivd 
viva 
viva 
viva 
viva 
viva 
wid 
viva 
vivo 
vivo 
vivd 
vivo 
vivo 
vivo 
vavd 
vivd 
vid 
viva 
vivd 
viva 
viva 
viva 
vivd 
viva 
vivd 
viva 
viva 
vivd 
viva 
vivd 
viva 


a | aah lates 


t (ore) orL 0°0 00S00S OS100S 720% biSy0S yuvad 
‘ 0°0 ork 0°0 00S00S SE!100S Z0z 01S00S uyvad 
‘ 0°0 o°k 0-0 00S00S O2100S Z0z% 60S500S uvad 
‘ 070 orn o°0 00S00S S0100S z02 B0S00S uveo 
iy 0°0 ort (ae) 00S00S 06000S 702 20500S yuvad 
t 0°0 ov4 0-0 00S00S SZ000S Zzoz 90S500S uvad 
b 0°0 ovt 0°0 00S00S 090005 Z0z% G0S00S yuvao 
‘ 0°70 ort ome) 00S00S SProsos zo0z vosoos uvad 
‘ 0°0 ovr 0°0 00S00S £00065 Zor €0S500S uvad 
ry 0°O on 0°0 00S00S St000S Z0z Z0S00S uvad 
ry 070 o's 0°0 00S00S o00000S 102 40S00S yvad 
$ 
SINIOd Y31N39 TVIIHLAW OL 3DYZANOD SGN3 MOBI3 L¥ S3xOdS $ 
$ 
wing NI93¢8 
33S 108 ONIOIAIOSNS STVAUZINI-Z 4O “ON SI & 
23S dOlL ONIGIAIOGNS SIVANZLNI-Z 3O “ON SI Pf 343HM 
wiser te assy 
HdI+(1«01+06S)*0001=01 (wol.0g) 
baP ices eo tg osy 
Hd I+( (1-0) *04-00S)*0001=G1 (OL) 
- (93S JHOITWH1S) 
IHd 430 S3N1VA IWGON G3QONNON 3yY HdI 
V13HL JO S3NTVA TVYGON Q3QNNOY 3¥V HLI 3Y3HM 
HdI+(HLI+00S)+0001=Q01 (m0873)°** S3NT1VA 1VvOON 


“h'd WON JONVISIG=Z (woiids) 


°I°d WOHS ZINVLSIG-=Z (dol) 
IHd=1 


3dlau=yu (IIS AHOTVELS)**"Sw3aLSAS ayo00d 


(THA INIS=3dTdt=Z 
(rata 'o'd LY O6°O=V13H1) ‘VLSHI=L 


(IHd )SO9*3d1ay¥-ONZQYN=a (MOG73)°**SW3LSAS GYxOOD 


(ogt‘O=IHd) WNDYID-37YH 3dId LV G3SOGWI**’ 


SNOT193S iHOSIveLS AB AaNIOPOY 
(O°St=QN2g4) mMOB13 33Y93IG 06 


30 SISISNOD (€91°S=3dI 


du) 3d 


lg'** 


AULIWW 


3yN1ING 


AS 


is 


AAAAAAARAAAARAAARAAAARAARAARAAAAAHRABAA 


"0997e° "0° °C’ O=03S8NS 
B1I-NI o99ZE€ JO IN3NOW TYNOISUOL 23 


VLATLENs 
be woosns 


“O0’o997f'°0* > O=03SENS 
@I-NI 099ZE JO INSWOW 3NW1d-40-1NO=3ITLILENS 


€2 woDEns 


“O° 'o' 0gaze-* °0=C4aSENS 
GI-NYI O997E 40 LN3WOW SNYIdG-NI=3T1 51 8NS 
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zon viva 0000%p» o0000r INYOND 


194s Wiwa S910ly o0000% 1NYDND 
oot =e viva OSLOly 000007 LNYDND 
6St ss wiva SE10lp 00000P LINYOND 
ests viva OZ710tp 000008 LNYOND 
2Sb—s Wiva SI10le o00000P LNYDND 
9S1 viva O600Ih 000007 INYOND 
ssi viva SLO0lp 00000¥ LNYOND 
vSt viva 09001r 00000r LNYOND 
est viva Sro0tp o0000P 1NYOND 
zwStswiva O£00!e 000007 INUSND 
Ske waa S1001b 00000% INYOND 
ost = viva 00001% 00000% LNYOND 
6bi = s viva S9100~ 000007 ANYOND 
epi = viva 0S100p oc000¥ SNHOND 
dys = wiva SE100p 00000¥ INYOND 
Ort = WG OZ100% 00000% LNYOND 
Spi Wava S0100y ocOc00Y INYOND 
tri viva O6000e 000054 LNYOND 
epi = yiva S!000% o0coor LNYOND 
tpi =e waa 09000» 00Cd0¥ LNYOND 
tyk = Wawa Sr000p 00000¥ LNYOND 
Ovi = viva C£000% 00C00r7 ANHYOND 
6€l viva $1000b. o00Cor INYDND 2 
set viva S9L06E O0000KE LNSOND 
zet = viva OSIOGE O0CDGE ANYOND 
9€b Viva SELOGE 0000HE LNYOND 
sti viva OZ106£ O0006€ LNYDND 
ves wiva GOLO6E O000KE INYDND 
e€€b viva OS006E O00006E 1NYOND 
zest = viva GLOOSE O00C06E LNYOND 
1€4 = Wiva OSOO0S6E 00006€ INYSND 
O£€i wiva : Srodse 00006E LNYOND 
6zi—s wiva OCO0S6E OOCOGE LNYOND 
ezt ss wiva GIOO6E O0006E INYOND 
£ts Ss waa $ 
9z1 viva SNOIL93NI0 HiO0S (93S iHOIvYLS) $ 
sts = vivo AINO NOILO3SNIO vVi3HL (M9813)°°°ADNSNYONOD 1N3W373 $ 
ets =e wiv $ 
€zi—s viva t 0°0 ob 0°0 00S06S oOBl06S 102 €1S06S uvao 
zis ava t 0°0 O'4 0°0 00S06S S$9106S 202 Z1S06S uveo 
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